Draft version April 21, 2010 

Preprint typeset using L£-T^X style cmulatcapj v. 6/22/04 



ON USING THE CMD MORPHOLOGY OF M67 TO TEST SOLAR ABUNDANCES 
Z. Magic and A. Serenelli 1 and A. Weiss 

Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany 

AND 

B. Chaboyer 
Dartmouth College, Hanover, NH, USA 
Draft version April 21, 2010 

ABSTRACT 

The open cluster M67 has solar metallicity and an age of about 4 Gyr. The turn-off mass is close to 
the minimum mass for which solar metallicity stars develop a convective core during main sequence 
evolution as a result of the development of hydrogen burning through the CNO-cycle. The morphology 
of the color-magnitude-diagram (CMD) of M67 around the turn-off shows a clear hook-like feature, 
direct sign that stars close to the turn-off have convective cores. VandenBerg et al. investigated the 
possibility of using the morphology of the M67 turn-off to put constraints on the solar metallicity, 
particularly CNO elements, for which solar abundances have been revised downwards by more than 
30% over the last few years. Here, we extend their work filling in the gaps in their analysis. To this 
aim, we compute isochrones appropriate for M67 using new (low metallicity) and old (high metallicity) 
solar abundances and study whether the characteristic turn-off in the CMD of M67 can be reproduced 
or not. We also study the importance of other constitutive physics on determining the presence of 
such a hook, particularly element diffusion, overshooting and nuclear reaction rates. We find that 
using the new solar abundance determinations, with low CNO abundances, makes it more difficult to 
reproduce the characteristic CMD of M67. This result is in agreement with results by VandenBerg 
et al. However, changes in the constitutive physics of the models, particularly overshooting, can 
influence and alter this result to the extent that isochrones constructed with models using low CNO 
solar abundances can also reproduce the turn-off morphology in M67. We conclude that only if all 
factors affecting the turn-off morphology are completely under control (and this is not the case), M67 
could be used to put constraints on solar abundances. 

Subject headings: Sun: abundances — open clusters and associations: individual (M67) 



1. INTRODUCTION 

The latest revisi on of the solar abundances by 
lAsplund et all (120051 AGS05) resulted in a drastic reduc- 
tion of the abundances of carbon, nitrogen, and oxygen, 
and thus in the total metallicity of th e Sun. While Z/ X 
in the previous sta ndard references bylGrevesse' fe Noels) 
(fl99l GN93) and lGrevesse fc Sauvall (fl998L GS981 was 
0.0245 respectively 0.0230. lAs plund et all (|2005| i de- 
termined a value of 0.0165. The rece ntly published 
compl ete re-analysis of all elements by Asplun d et al.l 
(2009, AGS09) corrected this number slightly upwards 
to 0.0183. It is interesting to realize that Z/X has con- 
stantl y gone down with time since lAnders fe Grevessel 
(|1989h . who gave a value of 0.0275. 

While the analysis by lAsplund et all (|2005| ) and in 
particular by lAsplund et al.l (|2009l ) is undoubtedly the 
most complete, coherent and technically advanced one, 
the results are far from being accepted. The reason is 
that standard solar models (SSM), calculated for the 
new metal abundances, agree much less with well es- 
tablished and accurately determined results from helio- 
seismology: the sound speed and density profile, the 
depth of the conve c tive en velope, and its helium con- 
tent. iBahcall et al.l (|2005al ) investigated the new "solar 
model problem" , which has been confirmed in many sub- 
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sequent, independent publications. lAntia fe Basul (2006) 
further demonstrated that helioseismic data require the 
solar envelope to have a co mposition close to th e old 
photospheric abundances, and lChaplin et all (|2007t ) used 
low-degree p-modes penetrating to the solar core to ar- 
rive at a similar conclusion. 

To cure the "trouble in paradise" (|Asplund et al.l 
2005) attempts were undertaken to confirm the previ- 
ous, higher abundances of intermediate elements, in par- 
ticular of oxygen (e.g. iCenteno fe Socas-Navarrol 120081 : 
Caffa uet al.ll20 08). On the other hand, the lower abun- 
dances were recovered, too ([Majorca et ail 120091 for ni- 
trogen; iSocas-Navarro fe Nortonll2007l for oxygen). 

On the side of the solar model community, possi- 
ble or necessary changes to the constitutional physics 
were discussed to restore the prev ious ex c ellent agree- 
ment with helioseismology (see Guzik| 120081 for a 
summary). In particular, an increase of opacities 
has been considered the most promising approach 
(iBasu fc Antiali 2004: Bahcall et al. 2005b; lAntiafe Basul 
120061: iChristensen-Dalsgaard et all 120091 ) . But updated 
Rosseland opacities for the solar interior a re not enough 
to sol ve the solar abundance problem (jBahcall et al.l 
l2005cf K Also, a postulated upward correction to t he neon 
abun dance has not found observational support (|Young| 
l2005h . 

It is therefore reasonable to shift attention also to other 
stellar objects, since any revision of the solar metallicity 
yardstick leads to a corresponding change in the metal- 
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licity of all stars. This is in particular true for stars with 
[Fc/H] as 0, since their metallicity is often determined in 
a strictly diff e rentia l and therefore quite accurate way. 
lAlecian et all (|2007l ) found that with the lower AGS05 
solar composition they could reproduce the properties of 
the pre-main sequence binary system RS Cha, while it 
was impossible to fulfill all observational constraints with 
the higher GN93 abundances. 

Another obvi ous test was performed by 

IVandenBerg et al.l (|2007l VG07), using the CMD 
morphology of the open cluster M67, which has a 
meta llicity of [Fe/H] = 0.00 ± 0.03 (|Randich et al.1 
20061) and a turn-off (TO) age around 4 Gyr (e.g. 



Pietrinferni et all [20041 . At this age the turn-off mass 



is slightly higher than 1 Mq and close to the critical 
mass at which a convective core on the main-sequence 
appears due to the dominance of the CN O-cycle over 
the pp- chain dKippenhahn fc Weigertlll990f). The CMD 
of M67 ([Montgomery et al.lll993l : lSandquistll2004[ ) shows 
a clear hook-like structure at the turn-off, indicative 
of a convective core for the turn-off stars. Since the 
efficiency of the CNO-cycle is directly proportional 
to the abundance of these elements, the reduced solar 
abundances of AGS05 could possibly result in a radiative 
core for stars at the TO, although in the fitting process 
the change in luminosity (at given mass, luminosity is 
lower for the AGS05 mixture) and effective temperature 
complicates simple predictions about the TO-mass 
and core structure. Indeed, VG07 found that the 
TO-hook disappeared for the new abundances, thereby 
strengthening the case for the GS98 metallicity scale. 

However, VG07 did not fail to mention some caveats: 
Their models were calculated without atomic diffusion 
taken into account, a physical process that at least for 
solar models - irrespective of the metallicity scale used 
- is essential for the best possible agreement with seis- 
mic inferences. They correctly pointed out that diffu- 
sion helps to support a convective core a t a stellar mass 
lower than in models ignoring this effect (Mi chaud et al.l 
2004), mainly due to the increase of CNO abundances 
in the stellar core by gravitational settling. The different 
ages of M67 and the Sun, but identical present-day abun- 
dances could also imply different initial abundances, an 
effect that should be taken into account in careful and 
precise studies. Apart from these points, one should keep 
in mind that with the occurrence of a convective core 
the question about the amount of overshooting arises. 
Its treatment, and in particular the value of any free pa- 
rameter in its practical implementation, cannot simply 
be taken for granted, as usually such "calibration" was 
based on stellar models resting on the higher solar metal- 
licity scale. Finally, the occurrence of a convective core 
may also depend on other aspects and therefore any con- 
clusion concerning the metallicity is valid only under the 
assumption that all other influences are under control. 

Given the result by VG07 and the open questions 
raised above, we revisited M67, trying to complete and 
extend the pioneering study of VG07. In § [5] we will in- 
troduce briefly the stellar evolution code we have used for 
most of the calculations in the present work, and those 
aspects of the input physics, besides the solar metallic- 
ity, we will also investigate. In § [3] we demonstrate that 
we can recover the results of VG07, a fact that is not 
irrelevant in the light of the following section. We then 



present our preferred model, showing that in this case the 
CMD morphology of M67 can be reproduced for both so- 
lar abundance scales. In § @] we investigate, how other 
approaches and variations of the input physics may also 
influence the quality of isochrone fits, and also possible 
systematic differences arising from the use of different 
stellar evolution codes. Our conclusions will follow in 

§E 

2. STELLAR MODELS AND CRUCIAL DATA 
2.1. GARSTEC 

Most of the stellar model calculations presented in this 
pape r were done with the GARSTEC (jWeiss fc Schlattll 
12008? ) program. Here, we briefly present aspects or mod- 
ifications of the code, which are relevant for this work. 
The standard assumptions and settings can be found in 
the quoted reference. 

Convective overshooting is treated as a diffusive 
process acco r ding to the approach described by 
iFrevtag et al.l (|1996f). The exac t implementation can 
be found in iHerwig et al.l (|1997| l Based on hydrody- 
namical simulations, the diffusive constant therein (D os ) 
is assumed to decay exponentially beyond the formal 
Schwarzschild-border as 

— 2z 1 
D os (z) = Do exp -j^-, D = -v ■ H P . (1) 
J tip 6 

This is equivalent to exponentially declining velocities 
of convective elements. In (TTJ) z is the radial distance 
from the Schwarzschild border and Hp the pressure scale 
height taken there. The constant Dg sets the scale of 
diffusive speed, and depends on the convective velocity 
i>o inside of the convective border. / is a free parameter 
defining the scale of the overshooting. It is known that 
fitting the CMD of young open clusters usually leads to 
an overshooting region extending for about 0.2 Hp in 
the classical local prescription. We have computed some 
test models for stars between 2 and 6 Mq and found that 
this overshooting value is well reproduced by the diffusive 
approach, during the hydrogen core burning phase, with 
a value of about / = . 018, a confirmation of previous 
results by IHerwig et al.l (|1997f ). 

For small convective cores, the amount of overshoot- 
ing has to be limited. In the local description this is 
done, for example, by a multiplicative factor, which in- 
creases linearly with m ass in the range of 1 to 2 Mq (see 
Pietrinferni et al. 2004, for an example). In our code this 
is achieved by a geometrical cutoff, where Hp in Equa- 
tion [T] is replaced by 



Hp := Hp ■ min 



fAR c 



V H f 



(2) 



where Ai? cz is the thickness of the convective zone. In 
this way, the geometric cutoff has the advantage that 
the overshooting region is always limited to a fraction of 
the extension of the convective region. It is our standard 
li miting procedur e for oversho oting. A similar cut is used 
in lVentura et~aT1 ([2005L fl998l ) . 

An alternative approach, namely the use of a ramp 
function for the / parameter, is also applied for some of 
our calculations. For masses lower or equal than 1.1 Mq 
we use / = (i.e. no overshooting), whereas for masses 
equal to or higher than 2.0 Mq we set f — 0.018. In the 
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intermediate mass range the efficiency of overshooting 
varies linearly with mass 

/ = (0.13M/M o - 0.098)/9. (3) 

The small discontinuity in the overshooting efficiency is 
of no consequence for our study. 1.1 Mq models with 
/ = 0.005, as would result from application of the above 
equation to this mass, lead, given the low overshooting ef- 
ficiency, to the same evolutionary tracks as models with- 
out overshooting. This prescription for the overshooting 
efficiency leads to evolutionary tra cks that reproduce well 
those bv lPietrinferni et all ([20041 ) . 

Atomic diffusion is treated within the same diffusive 
numerical scheme. In this paper, hydrogen, helium and 
heavier elements (including iron) are diffusing. Radiative 
levitation is not taken into account. 

To agree with the physical assumptions in VG07, the 
default set of nuclear reaction rates are those of the 
NACRE collaboration (jAngulo et al.lll999D . To inves- 
tigate their influence, for individual reactions alternative 
rates were used. For the 14 N(p , 7) 15 Q reaction, we also 
used the rate from lMarta et al.l ((20081 ) . the newest result 
from the LUNA collaboration, which is lower by about a 
factor of two with respect to the NACRE rate at the rele- 
vant temperatures. This has a considerable effect on the 
TO morphology. Another reaction that turned out to be 
of surprisingly strong influence is the 17 0(p, a) 14 N reac- 
tion. We test ed this by employing either the recent mea- 
surements bv lMoazen et al.l (|2007l) , which for T > 4 x 10 8 
K is very sim ilar to the NACRE recomm endation, or the 
one given bv lCaughlan fe Fowler! (|1988| ). A detailed dis- 
cussion on the reaction rates can be found in § 14.21 

For all solar composition choices (AGS05, GS98) con- 
sistent Rosseland mean opacity tabl es were prepared 
follow ing the procedure described in IWeiss fe Schlattll 
((20081) . 

The final step to compare isochrones with observed 
CMDs i s the transformation t o color s. For this we used 
that by IVandenBerg fe Cleml (|2003l VC03) which, to- 
gether with a choice for the distance modulus and red- 
dening of M67, results in satisfying CMD fits on the main 
sequence and subgiant branch. Alte r native ly, we used 
the transformations by ICassisi et al.l ((20041 ) for testing 
purposes. 

The initial stellar parameters (Yin, Z m , mixing length 
parameter) are obtained from solar model calibrations. 
This will be discussed in more detail in § 13.21 With these, 
we computed the evolution from the zero-age main- 
sequence (ZAMS) to the tip of the Red Giant Branch 
(RGB) for mass values from 0.6 to 1.5 Mq in steps 
of 0.1 Mq. We recall here that, as in VG07, we are 
not interested in considering the RGB of M67, where 
deficiencies in our treatment of the outer layers of the 
star could have an impact on the location of the mod- 
els in the HRD. We confirmed by tests starting from 
the pre-main sequence that for stars with mass below 
the critical mass for the occurrence of a convective core, 
Af ccc , the transient convective core at the end of the 
pre-main sequence phase also appears here. This con- 
vective core is the result of the short phase of CN- 
conversion, but may be sustained if convective overshoot- 
ing is included. Finally, for constructing isochrones, the 
tracks are normalized to the so-called equivalent points 
(jBergbusch fe Vandenbe7il992tlPietrinferni et al.ll2004l) 



and the interpolation to an isochrone is done between the 
normalized tracks. After an isochrone age was fixed, we 
calculated an additional model with the TO-mass and 
recalculated the isochrone to make sure that the TO- 
morphology does not depend on the isochrone interpola- 
tion scheme. 

2.2. Data for M67 

To compare our models with M 67 we used the photo- 
metric data by ISandquist] ((20041) . These are more accu - 
rate than the older data from I Montgomery et al.l (|1993h . 
which were used by VG07, and contain bona fide sin- 
gle stars only. Thus, this CMD has a narrower main- 
sequence band. However, we also tested some of ou r 
isochrones with the data bv iMontgomerv et al.l (|1993h . 
We used both {B - V) and (V - I) colors. 

The dereddened distance modulus to M67 is (m — 
M) v = 9.70 a ccording to VG07, and 9.72 ± 0.05 
(|Sandauistll2004D based on subdwarf fitting t o the lower 
mam-sequence following IPercival et al.l ((20031) . The red- 
deni ng is E(B - V) = 0.038 (VG07) in good agreement 
with lSaraiedini et al.l (|1999ft . who gave 0.0 4 ±0.02. Sim- 
ilar va lues have been recently obtained bv lTwarog et al.l 
(2009). For the metallicity we assumed a value of 
[Fe/H] — 0.0, which is well withi n spectroscopica lly de- 
termined errors, for example by iGrattonl ((2000), who 
gave [Fe/H] = 0.02 ± 0.06. 

2.3. Model composition 

To be consistent with VG07, we used for the old respec- 
tively new solar abundances those by GS98 and AGS05, 
although the former are simply an update of GN93 tak- 
ing into account additional literature values, and AGS09 
would be the most recent and complete re-analysis. Since 
AGS09 abundances are slightly higher than AGS05, our 
choice is testing the more extreme case. The initial abun- 
dances of the stellar models for the M67 isochrones are 
taken to be the same as those resulting from solar model 
calibrations. 

3. STANDARD FITTING PROCEDURES 
3.1. Recovering VG07 

The idea of VG07 to use M67 for testing the effect of 
the new solar abundances rests on the fact that the turn- 
off mass of this open cluster is very close to the critical 
mass for the onset of core convection , M ccc , or "tran- 
sition mass" ((VandenBerg et aT1l2007l ). This transition 
mass depends on the aver age exponent of the energy gen- 
eration rate in the core fKippcnha hn fc Weigertlll990L 
Chap. 22), that increases with the contribution from 
the CNO-cycle, for which ecNO ~ T 15 , in contrast to 
e pp ~ T 5 for the pp-chains at the relevant temperatures. 
The importance of the CNO-cycle depends on the nu- 
clear reaction rates (see § 14. 2|) as well as on the amount of 
"catalysts" , i.e. the sum of the CNO-abundances. These 
are, incidentally, the elements with the largest reduction 
in their abundance according to AGS05. Therefore, the 
morphological change in the CMD, i.e. displaying the 
characteristic hook at the TO, or a gap in star density, 
can be used to determine whether the TO-mass is below 
or above M ccc . 

However, the sensitivity of the M67 TO morphology 
implies that also other effects may influence it, apart 
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Fig. 1. — Isochrone fit to M67 for the two solar mixtures GS98 
and AGS05 using a solar model calibration as in VG07 (see Tablc[l} 
and NACRE nuclear reaction rates. Isochrone ages are indicated; 
distance modulus and reddening agree with VG07 



from the metallicity. The nuclear reaction rates were al- 
ready mentioned. Other possible aspects could be the 
amount of overshooting, atomic diffusion, and the pre- 
main sequence history. Last, but not least, technical de- 
tails of the stellar evolution codes may play a role. It is 
therefore necessary to show that our results do not de- 
pend on the particular numerical code. As a first, and 
crucial step, we have successfully attempted to repro- 
duce the key result by VG07. This will also show where 
changes to their procedure are indicated. 

The first step concerns the solar model calibration. 
Here, as for the M67 model s, VG07 used NACRE nu clear 
reaction rates, the O PAL (jlglesias fc Rogers! 119961 ) and 
I Ferguson et al.l (|2005f ) opacities (as is done in our code), 
but ignored diffusion. This point is crucial, as we will 
show later on. Some amount of convective overshooting 
(see § I4.1[) was included for the M67 models, but is not 
relevant for the solar calibrat ion. For the atm ospheres 
they used the MARCS model (|Gustafssonll2003l ) . In this 
respect we differ, since we use standa rd Eddington grey 
atmos pheres. However, as shown by IVan dcnBcrg et all 
(2008), this has in our present case no significant influ- 
ence on the tracks on the main-sequence and subgiant 
branch. 

Table Q] contains in the first two rows the resulting so- 
lar model parameters by VG07 and in rows 3 and 4 the 
equivalent ones obtained with our code. Additionally, 
we also show the results when using the Dartmouth and 
LPCODE codes (§ [O]). Note that the mixing length 
parameters can never be compared in their absolute val- 
ues due to the different formulations of MLT. Overall, 
the agreement with VG07 for the initial abundances is 



TABLE 1 

Results of solar calibrations for different codes and 
solar abundances (cols. 1 and 2) and for cases with and 

WITHOUT DIFFUSION (COL. 3). Y in AND Z in ARE THE INITIAL 
HELIUM AND METAL ABUNDANCES. A DENOTES THE RELATIVE 
DIFFERENCE WITH RESPECT TO VG07 FOR THE CORRESPONDING 
ABUNDANCE. 





Code 


Mixture 


Diffusion 


" M L 1 


1 in 


A 


"in 


A 




VG07 


GS98 


no 


1.84 


0.26760 




0.01650 








AGS05 


no 


1.80 


0.25590 




0.01250 






GARSTEC 


GS98 


no 


1.59 


0.26109 


-2% 


0.01661 


1% 






AGS05 


no 


1.62 


0.25018 


-2% 


0.01217 


-3% 






GS98 


yes 


1.72 


0.26903 


1% 


0.01866 


13% 






AGS05 


yes 


1.75 


0.25939 


1% 


0.01380 


10% 




Dartmouth 


GS98 


no 


1.78 


0.26461 


-1% 


0.01659 


1% 






AGS05 


no 


1.77 


0.24991 


-2% 


0.01219 


-3% 






GS98 


yes 


1.94 


0.27419 


2% 


0.01889 


14% 






AGS05 


yes 


1.91 


0.26075 


2% 


0.01405 


12% 




LP 


GS98 


no 


1.71 


0.26466 


-1% 


0.01648 


-0.1% 






AGS05 


no 


1.66 


0.25026 


-2% 


0.01217 


-3% 






GS98 


yes 


1.85 


0.27053 


-1% 


0.01806 


10% 






AGS05 


yes 


1.79 


0.25728 


1% 


0.01340 


7% 



within 1-2%, with our codes returning systematically 
lower initial metallicities (by 3%) for the AGS05 mix- 
ture. Although a small effect, it additionally disfavors 
the appearance of a convective core in the M67 TO. 

Using these initial composition values we calculated 
stellar tracks and isochrones for M67, which we show 
in Figure [1] The numerical values for distance and red- 
dening are identical to VG07. The best-fitting isochrones 
have somewhat higher ages (by 0.3 Gyr). The TO mass is 
1.229 respectively 1.200 M Q for the old and new compo- 
sition. One recognizes the same basic result as in VG07: 
for the AGS05 mixture, the isochrone does not show the 
characteristic hook. However, it displays a slight incli- 
nation indicating the presence of a very small convective 
core at the TO mass. Consistent with this morpholog- 
ical sign, the TO mass is marginally higher than M ccc , 
which is 1.175 M , while in the case of the GS98 mix- 



ture it is clearly above M c 



1.139 M m . These values 



appear to be very similar to those found by VG07 (see 
their Table 1). This exercise demonstrates that we are 
able to reproduce correctly VG07 and that their result is 
independent of the stellar evolution code used. 

3.2. Inclusion of diffusion 

VG07 already emphasized that atomic diffusion (effec- 
tively sedimentation) could modify this result as it leads 
to an increase of metallicity in the core over time. As 
explained in the previous subsection, this will favor the 
occurrence of a convective core and reduce M ccc . How- 
ever, the solar model calibration should also include dif- 
fusion, as it was shown that only with this physical effect 
the best agreement with the seismic Sun can be achieved. 
Only in this case the solar parameters, in particular the 
initial composition is determined as accurately as possi- 
ble. 

The corresponding results are listed in Tabled] as well. 
Obviously, both helium and in particular metallicity are 
now much higher as compared to the previous case ig- 
noring diffusion in the calibration; the latter values are 
increased by 10-13% (similar values for the comparison 
codes). We are using the same solar initial abundances 
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for the M67 isochrones. This is not completely consistent 
as M67 displays now a solar metallicity, but due to its 
slightly younger age diffusion should have started with a 
somewhat lower metallicity than the Sun. However, this 
effect amounts to a few parts in 10 in Z only, such that 
we did not iterate further the initial composition. 




(B-V) 

Fig. 2. — Isochrone fit to M67 using an initial composition ob- 
tained from standard solar model calibration with diffusion. Tracks 
for M67 include diffusion and NACRE rates; no overshooting was 
considered. 

The final fits to M67 are shown in Figure HJ We still 
used the same distance and reddening values, but the 
isochrone ages are now lower by 0.3 Gyr compared to 
FigureQ] The agreement with the ages by VG07 is purely 
incidental. Instead, the age for the G S98 mixture of 
3.9 G yr should be compared to that by Michau d et al.l 
(2004) of 3.7 Gyr. Now, for both compositions the CMD 
morphology is reproduced, and AGS05 is no longer dis- 
favored. The fit itself is of the same quality as before. 

4. VARIATIONS 

In this section we investigate the importance of fur- 
ther physical and numerical aspects that influence the 
appearance and size of a convective core on the main se- 
quence. These are obviously the nuclear reaction rates of 
the CNO-cycle and the amount of overshooting. Over- 
shooting may lead to the persistence of the otherwise 
transient convective core that may appear towards the 
end of the pre-main sequence, but may be kept alive 
due to additional mixing of CNO-nuclei by overshooting. 
With this aspect we start our tests. 

4.1. Overshooting 

We first repeated the models of § 13. 1[ i.e. diffusion 
was completely ignored, but included overshooting ac- 
cording to our standard prescription outlined in § 12.11 




0.4 0.6 0.8 1.0 1.2 

(B-V) 

Fig. 3. — As Fig.[T] but here overshooting with a parameter of 
/ = 0.018 and geometric cutoff was used. 



Equations [T] and [2J The parameter / can be varied, but 
we ke ep it below the standard value of 0.02 (|Herwig et al.l 
1997). In this work we take / = 0.018 because with 
this value, in our alternative approach (Equation [3]), we 
are able to reproduce w ell the evolutionary tracks from 
iPietrinferni et al.l (|2004t l , as mentioned in § 12.11 In Fig- 
ure [3] the resulting CMD for / = 0.018 and the stan- 
dard cluster parameters is displayed. We recall that the 
amount of overshooting is reduced due to our geometrical 
restriction (equation [2J ; this reduction can be quite re- 
strictive for stars with lower masses, around 1.2 M@, and 
be crucial in this study. However, even if the fit itself is 
not as good as before, now both mixture cases have TO 
masses with a convective core being present. 

The situation seems to be typical for cases without dif- 
fusion i n both the solar calibrat i on an d the cluster mod- 
els. In IVandenBerg fc Stetson! (|2004h the authors con- 
cluded that a small amount of overshooting is needed 
to reproduce best the CMD of M67. Ov ershooting 
was tr eated following the integral criterion by iRoxburghl 
(1989) and the adaption of a varying amount of over- 
shooting for v ery small convective core s close to M ccc , 
determined by IVandenBerg et al.1 (|2006f ). We stress the 
fact that this procedure was calibrated using open clus- 
ters with solar abundances on the old GS98 scale. For 
complete consistency, this procedure should be repeated 
for the new AGS05 composition, and this may explain 
why the same small amount of overshooting used by 
VG07 did not lead to a persistent convective core in case 
of the new abundances. 

The case of § 13. 2\ i.e. the one with full consideration 
of diffusion, was also extended by including overshoot- 
ing in the tracks for M67. We refrain from showing the 
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Fig. 4. — CMD of M67 with isochrones resultin g from tracks com - 
puted with the 14 N(p,7) ls O reaction rate bv IMarta et al.l IMPOST ) 
instead of that in the NACRE-library. In both the solar calibration 
and these tracks diffusion and overshooting were ignored. 



Fig. 5. — As Figure [4] but diffusion was included in both the 
solar calibration and in the tracks for the M67 isochrones. Left 
panel: overshooting constrained by the geometric cutoff (Equa- 
tion [2j; right panel: overshooting constrained by a linear depen- 
dence on stellar mass (Equation [3j. 



result here, as for both mixtures the convective core, as 
expected, persisted, but the overall fit g ot worse. This co- 
incide s with the isochrone fit shown bv lPietrinferni et al.l 
(2004), where the CMD could be better reproduced if 
overshooting w a s com pletely ignored. It agrees also with 
iMichaud et al.l (|2004f ). who neither found any necessity 
for overshooting, and who, too, did include diffusion in 
their models. 

4.2. Nuclear reaction rates 

As mentioned above, the nuclear reaction rates of the 
CNO-cycle are equally important for the occurrence of a 
convective core as is the abundance of CNO-nuclci. So 
far we have presented models employing the NACRE- 
library reaction rates for these reactions. The bottle- 
neck reaction of the CNO-cycle, 14 N(p, 7) 15 0, that de- 
termines the overall cycle rate, has been measured at 
stellar e nergies in the laboratory by the LUNA c ollab- 
oration (|Formicola et al.1 [20041 : IMarta et al.l 12008). and 
was found to be lower by about 50% with respect to the 
NACRE rate. The consequences for globular cluster age 
determinations and for some sensitive phases of low- and 
int ermediate- mass star evolu tion has been inve stigated 
bv llmbriani et al.l (|2004| ) and I Weiss et~aT1 (|2005l) . 

The rate being lower, we expect that the transition 
mass to harboring a convective core increases. There- 
fore isochrones using the lower CNO-abundances of the 
AGS05 mixtures will be even less likely to show the 
TO hook. We repeated the case of § 13.11 (i.e. ignoring 
diffusion completely) with the updated and most likely 
more accurate rate bv IMarta et al.l ([2008). The resulting 
CMD for M67, using again our standard distance mod- 



ulus and reddening i s show n in Figure [4] In agreement 
with llmbriani et al.l ([2004D the isochrone age had to be 
increased (by 0.3 Gyr for both mixtures) and M ccc in- 
creased by 0.08 Mq in both cases. As a consequence 
also the GS98 case now lacks the characteristic hook. We 
add briefly that the inclusion of overshooting does not al- 
ter this result because the geometric cutoff is restrictive 
enough that the convective region formed at the end of 
the pre-main sequence phase can not be maintained dur- 
ing the main sequence evolution. The values for Af ccc are 
1.215 (GS98) and 1.258 (AGS05), those at the TO 1.202 
resp. 1.196 M Q . 

This illustrates the fundamental problem with such 
tests: Imagine we would assume that the GS98 solar 
composition is the correct one, but we want to test which 
reaction rate is to be preferred. From Figures. [T] and 0] 
we would clearly conclude that the older one is to be 
preferred! 

We now repeat the computation of § 13. 2\ i.e. the case 
with diffusion included, but with the newer and lower 
14 N(p, 7) 15 reaction rate. In this case we recognize 
that with the older solar composition a small convective 
core is present (TO-mass/AT ccc = 1.214/1.172), which is 
completely absent in the AGS05 case (1.201/1.241 M Q ). 
While the CMD fit is not very good even with the GS98 
isochrone, it can be improved by including overshooting. 
The results, including overshooting with the geometric 
cut, are shown on the left panel of Figure [5j As men- 
tioned in the previous section, the geometric cutoff limit 
for overshooting is very restrictive particularly for lower 
masses. Therefore we recompute the tracks with the al- 
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Fig. 6. — Rates for the nuclear reaction 17 0(p, a) 14 N relative to 
that of the NACRE library, in the interesting temperature regime 
for core hydrogen burning in typical M67 stars near the TO mass. 



tcrnativc limiting algorithm, the ramp function described 
by Equation[3] Changes in TO morphology for the GS98 
composition are only minor. On the other hand, for the 
AGS05 composition, we now get a similarly good CMD 
fit as with the old mixture (right panel of Figure [5]), as the 
convective core that appears in the pre-main sequence 
can be maintained during main sequence evolution. 

We found a second reaction that is also influencing 
the efficiency of the CNO-double cycle, though indi- 
rectly. This is the 17 0(p, a) 14 N reaction that closes the 
CNO-II cycle and which is, together with the competing 
17 0(p, 7) 18 F the slowest one of that subcyc le. For this re- 
action new cross section measu rements by iMoazen et al.1 
(|2007[ ) and IChafa et al.l (|2007l ) are available. The re- 
sulting two new reaction rates are quite similar to each 
other and vary in the interesting temperature range of 
15 - 22 MK between 50% and 120% of the NACRE rate. 
Additionally, we have also considered for this reaction 
the rate given bv lCaughlan fc Fowler! ()1988l CF88). The 
comparison of all rates is shown in Figure [6] 

In Figure[7]we show the evolutionary track of a 1.2 M Q 
stellar model with GS98 composition for different choices 
of the two key nuclear reactions mentioned above. The 
track corresponding to our standard choice of nuclear 
rates, NACRE, is shown by the dotted line. Now, if the 
CF88 rate for the 17 0(p, a) 14 N reaction is used, which is 
less than 10% of the NACRE rate throughout the rele- 
vant temperature regime, the track for this crucial mass 
is strongly modified and no longer shows a sign of a con- 
vective core, as depi cted by the long-dashed line. If the 
IMoazen et all (2007) rate is used, the differences with 
respect to the track using NACRE rates are very small; 
for this reason, we do not show this track. However, it 
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In the figure, the legends denote the source for the reaction rates: 
left refers to the 14 N(p,7) ls O reaction, bottleneck of the whole 
CNO-cycle; right refers to the 17 0(p, o) 14 N reaction, bottleneck of 
the CNO-II cycle. 



is interesting to note that the effect that changing the 
17 0(p, a) 14 N rate has on the tracks depends on the rate 
used for the 14 N(p,7) 15 rate. Indeed, as stated before, 
using the LUNA rate for this reaction (and NACRE rates 
for all others) leads to the absence of a convective core 
in the evolution of a 1.2 M Q stellar model. This is shown 
in the short- dashed line in Figure [71 However, when in 
addition the Moaz e"n~et al.l (|2007[ ) rate is used, a small 
convective core is apparent in the track, as shown by the 
solid line in the same figure. 

We have identified the reason for this to be the follow- 
ing: the low 17 0(p, a) 14 N rate from CF88 slows down the 
CNO-II subcycle, reducing thereby the flux of 14 N back 
into the CNO-I cycle, which produces the overwhelming 
majority of the energy. While the loss of energy from 
CNO-II is rather unimportant, the drainage of available 
14 N by up to 90% from the CNO-I cycle is significant for 
its efficiency. One may consider this as a storage of nitro- 
gen in the form of useless 17 0. This makes the appear- 
ance of a convective core more difficult, and thus raises 
M ccc . Note that for the higher, older, 14 N(p,7) 15 rate 
from NACRE this effect is not important; it appears that 
in this case the branching is favoring the CNO-I cycle in 
any case, and the CNO-II cycle is always negligible. 

Since the importance of the 17 0(p, a) 14 N was rather 
surprising, and to make sure that our H-burning network 
included all necessary reactions, we did a test calculation 
with a variant of our code (Alves-Cruz 2009, private com- 
munication) that includes an extensive p-capture net- 
work up to silicon. In particular, the 17 0(p, 7) 18 F reac- 
tion, which is the branching reaction to the next higher 
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derived age is 4.5 Gyr, while for GS98 is 4.2 Gyr (for both codes). 



cycle, is included. We saw no difference in the evolu- 
tion of stars in this mass range up to post main-sequence 
phase. Our treatment of the CNO-cycles thus appears 
to be complete and correct. 

4.3. Comparisons with other codes 

We have seen that the behavior of the convective core 
for masses around the TO mass in M67 is very sensitive 
to the constitutive physics included in the stellar models. 
Given this sensitivity, it is important also to determine 
if calculation done with different stellar evolution codes 
lead to the same conclusions. Here we have used, in ad- 
dition to GARSTEC, results from two additional codes: 
LPCODE and the Dartmouth code. We summarize our 
findings below. 

LPCODE has originally been developed at La Plata 
Observatory. Exte nsive descripti o ns of the co de have 
been presented in lAlthaus et all (|2002l 120031 ). Here 
we only describe updates and changes done in the 
code that were specifically implemented for this work. 
The equation of state has been updated to the lat- 
est release by OPAL 2 . Radiative opacities from OPAL 
were calculated for both GS98 and AGS05 composi- 
tions and complemen t ed at low temperatures with those 
from iFerguson et "all ((2005). Conductive opacities, al- 
though not relevant for this work, are no w those from 
Potek hin and collaborators, presented in ICassisi et al.l 
(2007). For the nuclear reaction rates, NACRE is the 
stan dard choice but the LUNA rate for the 14 N(p,7) 15 
from IMarta et al.l (|2008f) has been adopted when neces- 

2 http://adg.llnl.gov/Research/ OPAL/EOS.2005/ 



sary. Overshooting is treated in a diffusive approach but 
the geometric cutoff (Equation [2} has been implemented 
to allow better comparison with GARSTEC results. 

We have computed a large set of models covering the 
ranges of mass, composition and input physics relevant 
to our work. In all cases, agreement between GARSTEC 
and LPCODE has been very good. Here, we choose to 
show in Figure [8] some evolutionary tracks that include 
the most complex input physics, as used for the CMD 
in in left panel of Figure [5] NAC RE rates upd a ted b y 
the 14 N(p,7) 15 reaction rate from IMarta et al.1 (|2008h . 
diffusion in both the solar calibration and in the tracks, 
and overshooting limited by our geometrical restriction 
(Equation [2]). Left and right panels show results for 
AGS05 and GS98 compositions, respectively. In both 
cases, tracks for 1.1, 1.2, and 1.3 M are shown, that 
bracket the TO mass of M67 and, consequently, deter- 
mine the morphology of the isochrones around the TO. 
By inspection of Figure [8] it can be readily seen that 
models with both codes show a very similar evolution 
along the HRD; differences are hardly noticeable. Re- 
sults are of a similar quality for other choices of con- 
stitutive physics and for both composition options. In 
relation to systematic uncertainties affecting the conclu- 
sions of our work, this result is particularly encouraging 
because GARSTEC and LPCODE have been developed 
completely independently, and do not share any numer- 
ical algorithm for solving the equations of stellar evolu- 
tion. 

The Dartmout h stellar evolution cod e is derived from 
the Yale code (|Guenther et al.l 119921 ). with modifica- 
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tions and updates descri bed i n IChabover et all ([ 2001), 
iBjork fc Chaboverl (|2006| ) and iDotter et al.l <|2008l ). For 
this project, the NACRE nuclear reaction rates were 
implemented, with the exception that the LUNA rate 
(jMarta et all 12008ft for the 14 N(p,7) 15 reaction was 
used. Convective core overshoot in the Dartmouth code 
is parametr ized as a multipl e of H p. Normally, the pre- 
scription of iDemarque et all ()2004j ) is used, in which the 
amount of core overshoot is small (0.05 Hp) for stars with 
small convective core masses, and gradually increased to 
0.20 Hp for stars with masses 0.2 M Q above the criti- 
cal mass for the turn-on of convection in the core. This 
prescription for convective core overshoot was found to 
yield good agreeme nt with open cluster co lor-magnitude 
diagrams, when the iGrevesse et all ()1996t ) solar mixture 
was used. This mixture has (Z/X) = 0.0244. For this 
project, a small constant amount of convective core over- 
shooting regardless of the stellar mass was used. The 
amount of convective core overshoot used in specific mod- 
els is identified in the figure captions, and was typically 
0.07 Hp. 

In the case of the comparison between results of 
GARSTEC and the Dartmouth code we choose to 
present CMD fits to M67 computed with both AGS05 
and GS98 compositions. Models for GARSTEC use over- 
shooting constrained by Equation |3l Results for the best 
fit isochrones for both codes are shown in Figure|9j where 
the left (right) panel shows results for the AGS05 (GS98) 
composition. As throughout this paper, our goal is to de- 
termine whether isochrones based on stellar models with 
one or the other composition can reproduce the occur- 
rence of the observed hook in M67. As can be seen in Fig- 
ure HI isochrones from both codes reproduce well the TO 
morphology for both compositions. In the case of AGS05, 
the age of for M67 obtained from both codes is 4.5 Gyr 
and for GS98 is 4.2 Gyr. Although the TO morpholo- 
gies are not identical for both codes, some differences are 
likely to be present since we have not attempted, even 
within the uncertainties given by the input physics and 
observational data, to obtain the best possible agreement 
between the codes and with the data. 

As a last comparison, we state that the results deter- 
mined in VG07 (see § I3.1j) could also be reproduced by 
both codes. 

From comparing results from GARSTEC with those 
from LPCODE and Dartmouth code, we conclude that, 
even if numerical schemes for solving stellar struture and 
evolution equations and implementation of physics are 
different in the different codes, our conclusions are ro- 
bust; they do not depend on the stellar code used. It 
implies that systematics between the codes are not an 
important source of uncertainty in our conclusions. 

5. DISCUSSION AND CONCLUSIONS 

The open cluster M67 has solar metallicity and its 
CMD shows a clear hook around the TO, evidence that 
stars populating the TO have convective cores. The mass 
of these stars is around 1.2 M Q , very close to the criti- 
cal mass value for which a convective core developes as a 
result of the efficient operation of the CNO-cycle. Based 
on this, VG07 have suggested that M67 can be used to 
test the new revisions of sol ar abundances presented b y 
Asplund and collaborators ()Asplund et al.l 120051 12009ft . 
since CNO elements, catalysts in the CNO-cycle, have 



suffered the largest reductions in their solar abundance 
values. Here, we have reconsidered the viability of using 
M67 as a benchmark for solar abundances, by studying in 
detail under which conditions stars with masses around 
1.2 M© develop (or not) a convective core towards the 
end of the main sequence and how this affects the mor- 
phology of isochrones best fitting M67 CMD. 

We have first focused on reproducing the results ob- 
tained by VG07. Using initial compositions for both 
AGS05 and GS98 obtained from calibrating solar mod- 
els without element diffusion, we have successfully re- 
covered their basic results: if stellar models include a 
small amount of overshooting but no element diffusion, 
then the hook in M67 is not reproduced by isochrones 
computed with the AGS05 solar composition, while it is 
present if the GS98 composition is used since, in this case, 
stars populating the TO have a well developed convec- 
tive core. This has been presented by VG07 as support 
for the GS98 solar composition; although VG07 pointed 
out to some caveats in their models that may influence 
this conclusion. 

We then included element diffusion in our models. On 
the one hand, solar models calibrated with diffusion yield 
a higher solar initial metallicity, as summarized in Ta- 
ble [1] Additionally, because gravitational settling is the 
dominant effect, metals tend to accumulate in the core. 
Inclusion of diffusion has then a double effect towards 
compensating lower metallicities. Our models with dif- 
fusion, but still using NACRE nuclear rates, lead to the 
formation of a convective core for TO stars in M67 even 
with the AGS05 composition and, consequently, the dis- 
tinctive hook observed in the CMD is also present in the 
isochrones, as shown in Figure [2] Under these assump- 
tions, we find no clear evidence favoring the older, higher, 
GS98 solar metal abundances over AGS05. 

In our study, we have also identified two nuclear reac- 
tions, 14 N(p, 7) 15 and 17 0(p, a) 14 N, that are critical in 
determining the precise value of M ccc . The first reaction 
is the bottleneck of the whole CNO-cycle, while the sec- 
ond one is the bottleneck of the CNO-II part. When the 
temperature is not high enough, the latter acts as a sink 
for 14 N nuclei, thereby slowing down the whole CNO- 
cycle. Let us first discuss results related to the first re- 
action. In addition to the NACRE rate, we also used for 
14 N(p, 7) 15 Q the la test measurement by the LUNA group 
(|Marta et al.l l2008f ) . which is about a factor of 2 lower 
than the NACRE rate. As can be expected, the LUNA 
rate shifts M ccc to larger values, by about 0.08 M . With 
this rate, even models with GS98 composition (not in- 
cluding diffusion) lack the characteristic hook observed 
in M67 (Figure |4]). This is not surprising in the light 
of results with AGS05 composition, since the effective 
operation of the CNO-cycle depends on the product of 
the abundance of catalysts and the 14 N(p, 7) 15 rate. 
The effect on the CNO-cycle of reducing the rate for this 
reaction by about a factor of 2, for stars with masses 
close to M ccc is very similar as to reducing the total 
number of catalysts by a similar amount. A somewhat 
more unexpected result relates to the second reaction, 
17 0(p, a) 14 N. The branching between the CNO-I and 
CNO-II cycles is almost independent of temperature for 
the temperature range that interests us in this work. 
However, if the temperature is not high enough for the 
CNO-II cycle to be fully active, 17 is created at the ex- 
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pense of 14 N but the feedback to the CNO-I cycle (the 
only one relevant in terms of energetics) is inefficient. 
As a consequence, the total number of catalysts is re- 
duced, with 14 N being stored in the form of 17 0. For 
this rate we hav e used, in addition t o the N ACRE rate, 
the older rate by Cauehlan fc Fowlerl ()1988l ) and the new 
measurement by Moazen et al.l (|2007f) . When we use the 
iCauehlan fc Fowled (11988) rate, that is about 10% of the 



NACRE rate in the relevant temperature range (see Fig- 
ure [S] for a comparison of the rates), we find that, even 
for GS98 and the NACRE rate for 14 N(p, 7 ) 15 there is 
no convective core in a 1.2 M Q models towards the end 
of the main sequence. This translates into the absence of 
a hook in the isochrones that would best fit M67 CMD. 
One can regard the last exerci se as merely of academic in- 
terest since the low rate from lCaughlan fc Fowlerl (|1988l ) 
seems now to be ruled out by experiments. However, the 
importance of this rate is still wor th being taken i nto ac - 
count with the newest results bv iMoazen et al.l (|2007h . 
even if th is rate agrees with N ACRE within 40%. By 
using the IMoazen et al.l (|2007| ) rate in the models, the 
1.2 M model with the LUNA rate for the 14 N(p, 7 ) 15 
reaction recovers a small convective core. These results 
are summarized in Figure [7J 

The development of convective cores in stars close to 
M ccc is sensitive to the detailed physical input in the 
models, as seen from the above results. Moreover, an 
appreciable amount of freedom in the modeling is still 
present because of the lack of a proper convection theory 
to account for, particularly, the amount of overshooting 
occurring in stellar cores. It is widely accepted, however, 
that an overshooting region of about 0.2 Hp is needed in 
stars with masses above ~ 2 M Q to explain, among oth- 
ers, the width of main sequences in stellar clusters. On 
the other hand, for stars below that mass, and particu- 
larly close to M ccc , it is also known that overshooting has 
to be limited to a smaller region. However, a detailed un- 
derstanding of how core overshooting depends of stellar 
parameters, namely the stellar mass, is not known. We 
have used here two different approaches to limit over- 
shooting, a geometric cutoff and an efficient factor that 
increases linearly with stellar mass, Equations [2] and [3] re- 
spectively, with the geometric cutoff being the standard 
choice. In relation to core overshooting, its importance 
is that it helps to reproduce the hook in M67 because it 
facilitates stars to sustain a convective core during the 
main sequence and the growth of the core towards the 
end of the main sequence for stars with masses around 
M ccc , critical for explaining the TO morphology in M67. 
On the other hand, overshooting cannot, by itself, cre- 
ate a convective core when there is none. This is the 
case, for example, for our models with AGS05 composi- 
tion and the LUNA rate for the 14 N+p reaction when the 
geometric cutoff is applied, as shown in the left panel of 
Figured] In this case, the convective core developed to- 
wards the end of the pre-main sequence evolution cannot 
be sustained by stars with masses close to the TO mass 
in M67, even with the help of overshooting. This may be 
the result of our geometric cutoff being too restrictive. 
When, instead, the linear ramp function of Equation [3] is 
used, isochrones based on AGS05 show again the charac- 
teristic hook of M67 (right panel in Figure [5]). With this 
prescription for overshooting, for stellar masses around 
M67 TO masses, the convective core survives during all 



main sequence evolution and naturally grows towards the 
end of it, due to increased CNO burning, leading to the 
occurrence of a hook in M67 CMD. Overshooting pro- 
vides a certain degree of freedom in stellar models, and 
here we show that this freedom renders isochrone fitting 
to M67 uncertain, in that depending on the way over- 
shooting is modeled, we are able to find, or not, a hook 
around the TO. It is important to note that both pre- 
scriptions for the limitation of overshooting for stellar 
masses around M ccc can be made consistent with cur- 
rent observational evidence. 

We have also tried to assess the possible dependence 
of our results on the numerics, that is the stellar evo- 
lution code. Therefore, we have used, in addition to 
GARSTEC, two other evolution codes: LPCODE and 
the Dartmouth code. None of these codes share a com- 
mon background in their development. For this work, 
when possible, we have tried to match the input physics, 
but some differences still remain. The most important, 
probably, is that in the Dartmouth code overshooting is 
accounted for as a fraction of the pressure scale height. 
For models appropriate for M67, this fraction has been 
taken as 0.07 Hp. Comparison of results between the 
three codes is highly satisfactory. Solar model calibra- 
tions, summarized in Table [1] yield very small differ- 
ences. A comparison between evolutionary tracks com- 
puted with LPCODE and GARSTEC is shown in Fig- 
ured] The agreement between results from the two codes 
is excellent. A comparison of isochrone fitting to the 
CMD of M67 with GARSTEC and the Dartmouth code 
is shown in Figure |H] for both compositions used in this 
work. Again, results are similar with both codes, which 
are able to reproduce the observed morphology around 
the TO. For each composition, both codes also yield the 
same age for the cluster, 4.5 and 4.2 Gyr for AGS05 
and GS98 compositions respectively. We conclude from 
these tests that systematics originating in the use of dif- 
ferent numerical schemes, and even some differences in 
the implementation of the input physics, do not affect 
our conclusions. 

Finally, we comment on two additional points that re- 
inforce our main conclusions given below. As stated in 
§ 12. 2\ the uncertainty in the metallicity of M67, accord- 
ing to iGrattonl (|2000h is about 0.06 dex, i.e. 15%. We 
have not played with the metallicity of our models and 
have always assumed a solar value (where solar means 
GS98 or AGS05). However, in view of the dependence of 
the precise value of M ccc to the details of the constitu- 
tive physics, variations in the metallicity of our models 
can in principle affect the occurrence or not of a convec- 
tive core for stars with M67 TO-mass. The last point 
is that in the new de termination of solar abundances by 
lAsplund et al.l ()2009D CNO elements have been revised 
upwards. The t otal n umber of CNO catalysts in the new 
lAsplund et al.l (|2009| ) is about 8% higher than in the pre- 
vious AGS05 solar abundance compilation. This increase 
would also contribute to facilitate the occurrence of a 
convective core at the TO-mass in M67. 

The main conclusions we draw from our study are as 
follow. The different solar compositions, namely AGS05 
and GS98 certainly have an impact on the predicted mor- 
phology of the CMD in M67. The occurrence of the 
TO-hook is more difficult to achieve for the low CNO 
abundances of AGS05. This is a confirmation of previ- 
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ous results by VG07. However, we have also found that 
other constitutive physics in the models, e.g. element dif- 
fusion, nuclear reactions, prescription of core overshoot- 
ing, also influence the stellar mass at which convective 
cores start to develop. This translates, in the case of 
M67, into isochrones that may have a hook around the 
TO even with the AGS05 composition or, on the con- 
trary, that show no hook even for the GS98 composition. 
There is a certain degeneracy in the constitutive physics 
underlying the presence of the hook in M67 that cannot 
be, at the present status of our knowledge, disentangled. 
The morphology of M67 could in principle be used to 
test solar abundances, but only under the strong assump- 
tion that all other factors affecting the TO morphology, 



more precisely the occurrence of a convective core in stars 
with masses around the TO-mass, are completely under 
control. Since we understand this is not the case, we 
conclude that M67 CMD morphology does not present 
a strong argument against low CNO abundances in the 
Sun. 
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